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I. INTRODUCTION 

Studies of the dynamical evolution of supernovae re- 
quire an accurate knowledge of the properties of nuclear 
matter over a large range of densities p, temperatures T, 
and proton concentrations P, 01 ■ The density values 
range from a few times the normal nuclear matter den- 
sity po to ^ 10~^/3o, the temperature may be as high as 
~ 20 MeV, and the proton concentration may even be 
close to zero. Based on the laboratory experiments, the 
properties of cold (T = 0) nuclear matter around the sat- 
uration density po with proton concentration limited by 
the stability valley of finite nuclei are more or less well 
known. Outside this narrow limit, the understanding of 
the properties of nuclear matter is fraught with many 
uncertainties. 

The density dependence of the nuclear incompressibil- 
ity and of the symmetry energy are among the key el- 
ements in the simulation dynamics of supernova explo- 
sion i. At supranormal densities, nuclear flow in ener- 
getic nuclear collisions aided by model-dependent calcu- 
lations helps to understand the nuclear incompressibility 
[Ij H, ■ Similarly, inference can be made on the density 
dependence of the symmetry energy from the comparison 
of theoretical predictions with experimental data on the 
differential flow of neutrons and protons and from the 
7r-/7r+, K°/K+ ratios, etc Disassembly of hot 

expanded nuclei offers one of the best tools to study the 
characteristics of the nuclear symmetry energy at subnor- 
mal densities [HI, ■ Experimental data related to 
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isotopic distributions [Tjl , isospin diffusion [ij, [T^ , 
and isoscaling [l^, constrain the density dependence 
of the symmetry energy in the subnormal region. There 
is considerable uncertainty, however, in all these extrac- 
tions. The importance of the nuclear symmetry energy 
can be gauged by the fact that higher symmetry energy, 
for example, leads to a lower electron (e~)-capture rate in 
the supernova collapse phase that may result in a strong 
explosive shock P, Tl8| . The isotopic abundance of rel- 
atively heavier elements in explosive nucleosynthesis is 
further directly correlated to the symmetry energy. 

Microscopic approaches based on realistic NN inter- 
actions, Brucckncr or variational schemes, or on ef- 
fective field theories show a large range of predictions 
on the density dependence of the nuclear incompress- 
ibility lii,!^, Jsinp, m, [13 and symmetry energy 
[3, HI, uS^K^^ These results are model depen- 
dent, which seems unavoidable at near-normal densities 
and above. At lower densities, however, based on the 
general analysis of the grand-canonical partition function 
for nuclear matter in the 5'- matrix framework [sol . [sij , it 
is possible to have predictions for various nuclear observ- 
ables which are practically model independent. For dilute 
nuclear matter, the system minimizes its total free en- 
ergy by forming clusters. The observables are expressed 
in terms of specific known properties of these cluster 
species like their binding energies, scattering phase shifts, 
etc. Using the virial expansion technique, Horowitz and 
Schwenk |32l | have evaluated the symmetry energy coeffi- 
cients of clusterized dilute nuclear matter where the clus- 
ter species were neutrons, protons, and a particles. In 
the iS-matrix approach, the calculations in nuclear mat- 
ter were extended with inclusion of all possible heavier 
clusters [l^l- The so-calculated results were found to be 
appreciably different from those obtained in Ref. [3^ . 
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in principle, may be infinite in size. In supernova mat- 
ter (charge- neutral due to the presence of electrons), the 
Coulombic term in the binding energies of the fragments 
severely constrains their size within the limits of the drip 
lines. The fragment composition is then likely to be 
altered, which would affect the properties of supernova 
matter in contrast to those in nuclear matter. The aim 
of the present paper is to investigate the symmetry and 
compression properties of supernova matter in the low- 
density regime (up to 0.02 fm~^) along with those for 
nuclear matter. Specifically, the density and tempera- 
ture dependence of the symmetry energy coefficients, the 
isovector potentials, and the symmetry incompressibili- 
ties are explored. 

The paper is organised as follows. The theoretical de- 
tails are presented in Sec. II. The results and discussions 
are contained in Sec. Ill, and the concluding remarks are 
given in Sec. IV. 

II. ELEMENTS OF THEORY 

The logical framework of the theory is set in the grand- 
canonical partition function of the interacting quantum 
system in the 5-matrix formalism of statistical mechan- 
ics as proposed by Dashen, Ma, and Bernstein [s^l- In 
subsection II. A, we present the key elements of the the- 
ory [sot applied to nuclear matter or supernova mat- 
ter. In subsection II. B, the expressions for a few rele- 
vant observables are given. In subsection II. C we present 
the methodology for calculating the single-particle and 
isovector potentials and the symmetry coefficients in 
some detail. 

A. The S-matrix framework 

The grand-canonical partition function for the two- 
component nuclear matter composed of neutrons and 
protons as the elementary species is written as 

Z = Tre^'^'-^^^p^p"^"^"^ , (1) 

where (3 is the inverse of the temperature T of the system, 
H the total Hamiltonian, iVp „ the number operators for 
protons and neutrons, and /ip.„ are the corresponding 
chemical potentials. The trace is taken over any complete 
set of states of all possible number of nucleons. Denoting 
the elementary fugacities by Cp = e'^'^" and Cn = e^^"^, 
the full trace can be decomposed as 

oo 

2= E iCpfiCnrTrz,Ne-''", (2) 

Z,N=0 

where Tiz,n is taken over states of Z protons and TV 
neutrons. For small Cp and C„, the quantity InZ can be 
expanded in a virial series 

\nZ = J2' DzMCpf iCnf ■ (3) 



Here the prime indicates that the term with Z = iV = is 
excluded. The knowledge of the virial coefficients Dz,n 
gives the partition function and thence the thermody- 
namic observables. 

In Ref. [ssj . it was shown that all the dynamical in- 
formation concerning the microscopic interaction in the 
grand potential of the system is contained in the partition 
fimction as two types of terms: 

In Z = In zj^J^^ + In Z^cat , (4) 

corresponding to contributions from stable single-particle 
states of clusters of different sizes (neutrons and protons 
included) formed in the infinite system and (multiparti- 
cle) scattering states, respectively. The superscript (0) 
indicates that the clusters behave like an ideal quantum 
gas. The first term can further be decomposed into con- 
tributions from ground states and excited states below 
nucleon emission threshold. Then, 

lnZ(l = lnZ(:!)+lnZ(?, (5) 

with 

Z,N •' ^ ' 

ln(lTC^.7ve-^(P^/2^™)). (6) 

The upper and lower signs in the above equation corre- 
spond to bosons and fermions, respectively, with go as the 
ground-state spin degeneracy. For mass numbers ^ < 8, 
the go values are taken from experiment. For heavier nu- 
clei, Qq is taken as 1 or 2 depending on whether the nuclei 
are bosonic or fermionic. The sum in Eq. extends over 
all possible fragment species that can be formed; Z, A^, 
and A refer to the proton, neutron, and mass number of 
a species, respectively. Here, p refers to the momentum 
of the fragment, m is the nucleon mass, and V is the 
volume of the system. The ejfective fugacity is given by 
Cz,N = e'^^^^'"+^^'"-', where the chemical potential of 
the fragment is fiz,N ~ Z^p + Nfi„ from the conditions 
of chemical equilibrium. The quantity Bz,n represents 
the binding energy of the fragment. 

In nuclear matter, only nuclear forces contribute to 
the binding energy as the Coulomb interaction is absent; 
then the fragment size can, in principle, be infinite. In 
supernova matter, the Coulomb interaction is operative 
in determining the binding energy and the size of the 
nuclei; the sum in Eq. (6) is then finite. The nuclei are 
embedded in a sea of electrons and the binding energy 
gets dressed up; the correction can be estimated using 
the Wigner-Seitz approximation (3^ . We work in natural 
units h = c = 1. The integral in Eq. (6) can be expanded 
in powers of the effective fugacity (^z,n as 
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The quantity = y^27r/ (AmT) is the thermal wave- 
length of a species of mass Am. 

A nucleus in a particular excited state is taken as a dis- 
tinctly different species and can be treated in the same 
footing as the ground state. The density of states is quite 
high in relatively heavy nuclei and increases nearly ex- 
ponentially with the square root of the excitation energy 
E* . Thus, the contribution of the excited states of a sin- 
gle nucleus is given as an integral over E* of the ideal 
gas terms weighted with the level density u!{A, E*): 




The expression for the level density is obtained from 
the Fermi gas model of nonintcracting nucleons in a nu- 
cleus im 

Lu(A,E)= . (9) 

The level density parameter a is taken as A/8 MeV"\ its 
empirical value. In Eq. (8), the prime indicates exclusion 
of light nuclei {A < 8) in the sum. The lower limit Eq 
is dictated by the location of the first excited state. We 
take it to be 2 MeV. The upper limit Eg is the particle 
emission threshold taken as 8 MeV. 

The scattering term in Eq. (4) can be formally written 
for the system under consideration as 

(10) 

where the sum is over all scattering channels, each having 
its chemical potential ^ and formed by taking any num- 
ber of particles from any of the stable species. The trace 
is over all the plane wave states for each of the chan- 
nels. 5" is the scattering operator and A is the boson 
symmetrization or fermion antisymmetrization operator. 
The subscript c, in diagrammatic language, refers to only 
the connected parts of the expression in parenthesis. 

To recast Eq. (10) explicitly in the context of nuclear 
or supernova matter, a set of channels with total proton 
number Zt, neutron number Nt, and mass number At will 
be considered. All other labels required to fix a channel 
are denoted by a. Obviously, the total mass and the 
chemical potential are independent of cr, depending only 
on Zf and Nt. The nonrelativistic energy in a channel is 
given by 

p2 

Ezt,Nt,cr = TTl Bzt,Nt,cr + i, (11) 

ZAtUi 

where Fcm is the total center of mass momentum, e is 

_ i_; _ i„ _ /^T\ J 3 7-1 1- _ 



sum of the individual binding energies of all the frag- 
ments in the channel. Integrating over Pqai, one then 
gets 

j-oc -pe / 8 \ 

Jo '^'^T^''Zum^\-^S-'{e)-S{e)j , 

fl2) 

the trace being now restricted to the channel {Zt, Nt,a). 
Examination of Eq. (12) shows that larger binding ener- 
gies are more important, because of the factor e'^^^' ™* " . 
Furthermore, two-particle channels are expected to be 
more dominant than the multiparticlc channels with the 
same Zt and Nt from binding energy considerations. The 
two-particle scattering channels are therefore only con- 
sidered. It becomes convenient to divide the channels 
into light ones, consisting of low-mass particles {A < 8, 
say) and heavy ones, containing at least one high-mass 
particle {A> 8), so that we write 

\nZscat = InZf^^f +lnZf„t. (13) 

The scattering of relatively heavier nuclei is known to 
be dominated by a multitude of narrow resonances near 
the continuum threshold. The S'-matrix elements are 
then approximated by resonances. Each of these reso- 
nances can be treated [s^ [13] like an ideal gas term. 
Then, In Z^^^ can be written in the form of In zi? , as- 
suming the resonance level densities to be the same as 
those of the excited states given by Eq. (9). The sum 
of the contributions from the excited and the resonance 
states can then be written as 

lnZ(? H-lnZf,, = T^E' r^E* u;{A,E*) x 

= ^E'>^(/lCz,A.±/2 If + •••). (15) 

The integration in Eq. (14) extends up to Er, the limit 
of resonance domination. The damping of the integral 
in Eq. (14) due to the presence of the Boltzmann factor 
limits the contributions to only those from low energies; 
we take Er= 12 MeV. The A-dependent entities 

fn{A) = dE*io{A,E*)e-^P''', (16) 

JEo 

with n = 1,2, decrease steadily with increas- 
ing n, so that the series converges quite fast. For 

the evaluation of InZ^^j (i-e., the contribution of 
j_ _ i„ ^ \ 1 j_j 1 _i i_ 
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NN, Nt,Nllc^,Na, and aa are considered, where N and 
t refer to the nuclcon and the triton, respectively. Then, 



hi Zf 



In +ln 



Each of the terms in Eq. (17) can be expanded in the 
respective virial coefficients. We consider the expansion 
up to the second-order coefficients which are written as 
energy integrals in terms of the relevant phase shifts. 

In summary, the grand partition function for the in- 
teracting nuclear system is given as 

InZ = InZ(J') + flnZ(o) + InZf,,) + In Zi^,„, . (18) 



Once the partition function is known, the chosen observ- 
ables can be calculated. Expressions for them in some 
detail are given in the following subsections. 



B. Equation of state 



The expression for In Z is given by 
InZ = F{ A[c„ + C, + ^-^a + ^-fd + Ib^pCnQ 



A? 



[Ct + 2CtibptCp + bntCn)] 



+ T3 [C/i + 2(h{bphCp + bnhCn)] 

+ i[C« + bc^c^C + 6a«Ca(Cn + Cp)]} + In Zfle.t • (19) 

In Eq. (19), the subscripts N, t, h, and a refer to the nu- 
cleon, triton, He^, and He'*, respectively. The coefficients 
b are the virial coefficients. In the limit of isospin symme- 
try, we take = bpp. Due to lack of pt scattering data, 
we assume bpt — bnh- These virial coefficients can be 
written in terms of experimentally known phase shifts. 
As an example, we write below explicit expressions for 
bnn and bnp'. 



bnn 



1 V2 

¥7^ + ^ 



/ dEditiE)e-P''^', (20) 



and 



with 



and 



bnp — b^np + bd : 



bd = 6V2 e^"/"^ 
V2 



r) f 



-I3E/2 



(21) 

(22) 
(23) 



In Eq. (21), the term corresponds to the non- 
resonance n~ p scattering contribution; the term bd cor- 



bound state of the deuteron with binding energy Bd- The 
energy E is measured in the laboratory frame. The ex- 
pression for the total phase shift is given as 



(17) S'^'^ = + = 0) + Sf;\l = 1)}.(24) 



LSJ 

The contributing partial waves are determined by the 
isospin / with the requirement of the antisymmetry on 
the total wave function of the NN system. 

The last term in Eq. (19) is the sum of the contribu- 
tions from the rest of the species {A > 4) and is given 

by 



lnZ«est = V'E;|(5o + c^{E*)e-'''/^dE*\.{25) 



As already stated, the second term in Eq. (25) con- 
tributes only for A > 8. The pressure can be evaluated 
from 



P = TlnZ/V. 



(26) 



The number density pi of the i-th fragment species is 
calculated from 



Pi = Ci 



d InZ 



(27) 



V,T 



The total neutron, proton and baryon density in the sys- 
tem can be obtained from 

i 

Pp = ^Z,p,, 

% 

p^ = Y,Mp^,. (28) 

i 

From the Gibbs-Duhem relation, the free energy den- 
sity is 



J" = -P + ^ PiP^ 



The entropy density S is calculated from 



(29) 



(30) 



which then yields the total energy density as 

Etot ^T + TS. (31) 
The detailed expression for the energy density is 



- E^*-^* - 0^{T3-(^nnCn + bnnCp + bnpCnCp) 



-j l^ j_; 
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+ -^Ct{bntCn + bptCp) + -JjChibrihCn + bphCp) 
1 

+ Tj[&«„Ca(Cn + Cp) + bac^Cl]} + T^iCpC^Cp 

+C(C' + Cp)} + -rrCtib'ntCn + b'ptCp) 
4T^ 

+ ^^Oi(^n/iC« + &p/iCp) 

+ TTRaCa+^a„Ca(C«+Cp)}, 



(32) 



where X^iG-ff denotes that the sum runs over the channels 
of heavy particles. In Eq. (32). the first term is identified 
with £cM, the kinetic energy density associated with the 
center of mass of the fragments. The second term refers 
to f *, the sum of the densities of the thermal and reso- 
nance excitation energies of the fragments, and the third 
term coming from the negative of the sum of the frag- 
ment binding energies Bi is denoted as £be- The rest is 
designated as Ef^, the sum of the contributions coming 
from the interactions between different pairs of light frag- 
ments. The primes on the virial coefficients denote their 
temperature derivatives. They are part of the entropy 
contributions. Thus, Eq. (32) can be rewritten as 



St. 



CM 



- £* + £be + Ej 



LL 



(33) 



C. The single- particle and isovector potentials and 
the symmetry coefficients 

The single-particle potentials for neutrons or protons 
in a nuclear medium are conventionally defined as [ssj 



(34) 



Here, Vo = [Vn -\- Vp) / 2 \s the isoscalar potential, Vi is 
the measure of the isovector potential, and r = ±i for 
neutrons or protons. The quantity X ^ iPn ^ Pp) I 
the asymmetry parameter, with as the total neutron 
or proton number densities in the system, and = p!^ + 
p^^ is the total nucleon density. The isovector potential 
V,sov is then 



Vn-Vp = ViX. 



(35) 



In the nucleonic medium, the single-particle potentials 
are obtained as 



Vr 



( dSi 
\dp?)pB^ 



(36) 



where £i is the interaction energy density. In homoge- 
neous nuclear matter, it is calculated from the effective 
interaction. In clusterizcd nuclear matter, we employ 
Eq. (36) for the definition of the effective single-nucleon 
potential. The calculation of Ej is, however, not straight- 



take counts of the contributions to the interaction en- 
ergy density. 

In Eq. (33), the first term is purely kinetic and the last 
term comes solely from the interactions. The binding and 
the excitation energies, however, are admixtures of both 
kinetic and interaction contributions. The binding term 
£be can be split as 



Ebe - - / 



^B,^-yp,{Bi 



K 



Bi). (37) 



In Eq. (37), the two terms are the kinetic and the inter- 
action parts of the binding energy in a fragment, respec- 
tively. The kinetic terms are estimated with the choice of 
the extended Hu'lthane wave function for the deuteron 
[sst and a Gaussian for t, He'^, and a particles [l^. The 
kinetic terms for these light nuclei are very small com- 
pared to the total binding term of the system, hence, 
even a change by 50% in the estimated values of the ki- 
netic terms in the light nuclei docs not change the values 
of the calculated observables noticeably. For heavier sys- 
tems, the Fermi gas approximation is employed. Since 
the binding energies are known experimentally, the in- 
teraction contribution Bl can then be obtained. Estima- 
tion of the interaction contribution £j from the excitation 
part £* is not easy. However, compared to the binding 
term ^ • PiBf , it is quite small. It is checked that changes 
in the single-particle and isovector potentials are insignif- 
icant with £j=0 or £j = £* . The calculations reported 
correspond to £j = ■ The interaction energy density 
then takes the form. 



(38) 



where we recall that £f^ refers to the contribution 
from the interactions between different pairs of light 
fragments. Using Eqs. (27), (36), and (38), the effec- 
tive single-particle and isovector potentials in clusterizcd 
matter can then be calculated. The single-particle poten- 
tials are called effective in the sense that they represent 
the average of the interaction of a single nucleon with 
all other free nucleons and bound nucleons in clusters 
of different sizes, the clusters being present in different 
proportions in the inhomogencous system. 

The total symmetry energy coefficient C*"* is defined 
through 



e{X) = e{X = 0) + Cf*X2 



(39) 



where e = £tot/ P^ is the energy per nucleon. The coeffi- 
cient C*°* can be split into interaction and kinetic terms 
as 



C*°* = Ci + C. 



K 



(40) 



corresponding to 
ei{X) 



ei{X = Q) + ClX^ 
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In nuclear matter, the isovector potential is related to 
the interaction energy per nucleon e/(X) as Visov = 
2{dei{X)/dX)pB. For homogeneous matter, e/(X) is lin- 
ear in X^, then 

V,sov = 4Cf X. (42) 

For clusterized matter, this relation is. however, only ap- 
proximate because ei{X) (as well as the energy per nu- 
cleon e{X) [3l|) is seen to be not fully linear in X^. 

The incompressibility K{X) of a nuclear system is cal- 
culated from 

KiX) ^ 9^, (43) 

where P, given by Eq. (26), is the pressure of the system 
with asymmetry X . It can be written as 

K{X) = K{X = 0) + KsX^, (44) 

where the coefficient A'g is the symmetry incompressibil- 
ity. 

III. RESULTS AND DISCUSSIONS 

We have calculated the equation of state (EOS), the 
single-particle and isovector potentials, symmetry energy 
coefficients and the symmetry incompressibility of dilute 
nuclear and supernova matter at different temperatures 
and densities. Results for nuclear matter are presented 
in subsection A and those for supernova matter are given 
in subsection B. Calculations have been restricted up to 
a baryon density = 0.02 fm"'^. At relatively high den- 
sity, the asymptotic wave function may not have a precise 
meaning and then expressions of the partition function in 
terms of ^-matrix elements may not be very meaningful. 

The virial coefficients related to NN, Na, and aa scat- 
terings along with their temperature derivatives are taken 
from Ref. [S^]- The same for Nt and Nh scatterings are 
obtained from Ref. In doing so, appropriate care 

has been taken for the slightly different choices of factors 
in the virial expansion of the partition function. 

A. Nuclear matter 

In Fig. 1, the pressure P is displayed as function of 
baryon density for symmetric {X~0.0) and asymmetric 
{X~0.3) nuclear matter. The upper panel corresponds to 
a temperature T= 4 MeV, the lower panel corresponds 
to T= 8 MeV. The results calculated in the S'-matrix 
approach (SM) are compared with those in a mean-field 
(MF) model. The MF calculations are performed with 
the SkM* interaction. At very low densities, P is nearly 
model independent {P ~ pT); at a little higher density, 
the difference between the two models is apparent from 
the figure. In the MF calculation, with isothermal com- 



beyond a certain density which increases with temper- 
ature. In the SM approach, with compression, the un- 
physical behavior does not arise because of many-body 
correlations (condensation). For symmetric matter, at 
low temperature, the pressure levels off at very low den- 
sities as shown by the full line in the upper panel sig- 
naling a behavior like a first-order phase transition. It 
points out a phase coexistence between light and heavy 
clusters. For asymmetric matter, the rise in pressure is 
linear in density at constant temperature for very dilute 
systems. There is a break from this linearity at a certain 
density (we refer to this density as the condensation den- 
sity); the pressure thereafter rises monotonically behav- 
ing like a second-order phase transition. At the higher 
temperature T= 8 MeV, the said transition occurs at 
a much higher density which is beyond the density we 
consider. The increase in free neutron multiplicity with 
density is mainly responsible for this monotonic rise. For 
both symmetric and asymmetric systems, the chemical 
equilibrium conditions coupled with the conservation of 
the baryon number and isospin governs this behavior. 

Since there is no Coulomb in the binding energy of the 
fragments formed in nuclear matter, the sum in Eq. (3) 
runs up to infinity in principle; in practice, one takes a 
finite sum for calculational facilitation. The calculations 
here have been performed with a maximum fragment 
mass Amax = WOO. The results are not very sensitive to 
further increase in the maximum mass (3l| . The binding 
energies of these nuclei are obtained from the liquid-drop 
type mass formula [4l| with Coulomb switched off. 

The effective single-particle potentials Vn and Vp for 
neutron and proton are shown in Fig. 2 as a function of 
density at temperatures T=A and 8 MeV for symmetric 
and asymmetric (X=0.3) nuclear matter. The isovec- 
tor potential Vn ^ Vp is also shown in the bottom panels 
for asymmetric matter at the same temperatures. Cal- 
culations have been performed in both the SM and MF 
approaches. In the mean-field model, the single-nucleon 
potentials decrease with increasing density with a mono- 
tonic increase of the isovector potential. They are nearly 
independent of temperature. In the S'-matrix approach, 
the single-nucleon potentials show a subtle behavior with 
density. For symmetric matter, Vn (—Vp), at lower tem- 
perature decreases sharply upto the condensation den- 
sity, beyond which it remains practically constant. Be- 
yond this density, the fragment composition of matter 
scales nearly with density and the constancy of the single- 
particle potential is a reflection of that. The value of this 
constant effective single-particle potential is ^ —34 MeV 
in contrast to that of '--^ —60 MeV for saturated uni- 
form nuclear matter. This is so because in calculating 
6£j/6p for uniform matter, one has to count the change 
in nucleon density, but in clusterized matter, the internal 
nucleonic density of fragments remains unaltered, only 
their number density changes. At lower temperature, 
for asymmetric nuclear matter, Vn goes through a sharp 
minimum around the condensation density. Beyond this 
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slowly; Vp on the other hand behaves more like that of 
symmetric nuclear matter, but it is deeper. The above 
qualitatively different behavior of Vn and Vp comes from 
the presence of free neutrons in asymmetric nuclear mat- 
ter. The isovector potential passes through a maximum 
and then decreases slowly with density. As in mean- field, 
the isovector potential is seen to be nearly proportional 
to the asymmetry X. The different nature of the single- 
particle potentials at higher temperature as seen for T~8 
MeV is a manifestation of the dilution of condensation 
effects. 

The symmetry energy coefficients, function of 

baryon density, are displayed in the left and right pan- 
els of Fig. 3 at temperatures T =4 and 8 MeV, re- 
spectively. In the MF approach, the total symmetry 
coefficient C*°* alongwith the kinetic and interaction 
components Cf and (shown as dot-dash lines) in- 
crease linearly with density and are seen to be practi- 
cally independent of temperature. The symmetry energy 
esym{X) = e{X) — e{X ~ 0) in the SM approach, as 
opposed to that in the MF model, is foimd to be anhar- 
monic in X , particularly at lower temperatures; this was 
already noted earlier [3l|. This induces an asymmetry 
dependence in the symmetry coefficients defined through 
Eqs. (39) and (41). To avoid this, we therefore, take the 
definition [sl] 

This applies to Cj, Cf^ and C*°* with appropriate choice 
of the energy components. The symmetry coefficients, so 
defined are shown as full lines in the figure. The symme- 
try coefficients thus obtained are seen to be very different 
from those obtained in the MF model. The magnitudes 
of the coefficients are much larger, the kinetic compo- 
nent is always negative and there is a marked dependence 
on temperature. The negative symmetry kinetic energy 
CfX"^ (= CKiX) -exiX = Qj) looks antiintuitive. For 
clusterized nuclear matter, it, however, can be under- 
stood from the fact that for symmetric matter, clusteri- 
zation is more favored leading to larger internal kinetic 
energy compared to asymmetric matter where there are 
more free neutrons and have lesser total internal kinetic 
energy. 

The anharmonicity of the symmetry energy esym{X) 
(calculated in the SM approach through Eq. (32)) in the 
asymmetry parameter is portrayed in Fig. 4. In the figure 
the evolution of Dse, a measure of this anharmonicity, 
is shown as a function of baryon density at different 
X at temperatures T =4 (full line) and 8 MeV (dashed 
line), respectively. The anharmonicity parameter Dse is 
defined as 

Dse ^ '^'^'''lZ:^T^''^ X 100 . (46) 

It represents the percentage deviation of esym{X) from 
C*°*X2 where Cf* is defined through Eq. (45). As ex- 

J 1 1 1 "j-l- J /^^_ _ 



other hand, it decreases with increasing temperature be- 
cause of the dissolution of clusters with heating. 

In Fig. 5, the incompressibility coefficient K of sym- 
metric (X=0.0) and asymmetric (X=0.3) nuclear matter 
are presented in the left panels as a function of density 
at T = 4 and 8 MeV. At very low density, the incom- 
pressibility in both the SM and MF approaches is nearly 
the same. Increasing density and subsequent condensa- 
tion renders the system more compressible in the SM 
approach. The somewhat different behavior of the in- 
compressibility for symmetric and asymmetric matter at 
lower temperature in the ^-matrix approach can be eas- 
ily understood from the different functional dependence 
of pressure with density for symmetric and asymmetric 
systems as shown in Fig. 1. In the MF approach, on 
isothermal compression, the system enters the unphysi- 
cal region beyond a certain density where the incompress- 
ibility becomes negative. This is not shown in the figure. 
From the right panels of this figure, it is seen that in the 
physical region, at a fixed density, the incompressibility 
K decreases almost linearly on isochoric cooling. In the 
SM approach, on isochoric cooling, there is a linear de- 
crease in K , however, a sudden fall in incompressibility 
at a certain temperature is noticed. This is related to 
the onset of condensation where the fragment multiplic- 
ity suddenly drops [4^ . 

The symmetry incompressibility defined in Eq. (44) 
is compared in the SM and MF models in Fig. 6 at tem- 
peratures T=4 and 8 MeV. As opposed to the MF model 
where Kg increases monotonically in the density region 
we explore, the same in the SM approach has a peaked 
structure which is prominent at lower temperature. Here, 
the symmetry incompressibility is seen to be weakly de- 
pendent on the asymmetry of the system while in the MF 
model it is asymmetry independent. The arrows in the 
figure indicate the density above which the system enters 
the unphysical region in the mean-field model. 

In Fig. 7, the symmetry incompressibility Kg is shown 
as a function of temperature at baryon densities 
0.001 and 0.01 fm~^. In contrast to the mean-field results 
where Ks is practically independent of temperature, in 
the S-matrix approach, is sharply peaked around the 
condensation region. The asymmetry dependence of Ks 
in this approach washes out with increasing temperature. 
The arrow in the lower panel of this figure marks the 
temperature below which the system is in the unphysical 
region in the MF model. 



B. Supernova matter 

Neutrons, protons and electrons are assumed to be the 
basic constituents of the supernova matter we are dealing 
with. The matter is also assumed to be in /3-equilibrium. 
For charge neutrality, the electron number density = 
Pp^ the proton number density. This /3-equilibrated mat- 
ter is treated both in the mean-field (MFB) and the S- 

J • /nAT■^»^ 1 mi^ _ J j_i ;„ _ 
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system depends on its baryon density and temperature. 
In supernova matter, in the SMB approach, the frag- 
ments have coulomb contribution in the binding energy; 
so, in the sum in Eq. (3), only nuclei within the drip 
lines are considered. For this purpose, all the isotopes 
(around 9000 in number) and their binding energies are 
taken from Ref. [i^ . The effect of electron environment 
is taken into account through the dressing-up of the frag- 
ment binding energies in the Wigner-Seitz approximation 
[3^ . [4^ . For a given electron density pe, the effective 
binding energy of the i-tli fragment species with Ai and 
Zi as its mass and atomic number is given by, 

B^ipe) = B,iO) + AB,, (47) 

where 
with 

- (-#)"'■ («) 

In the above equations, po is the saturation density of 

1/3 • 

normal nuclear matter, Roi= 1.16 A/ fm is the radius 
of the fragment of mass Ai and charge Zi. For proton, i?oi 
is taken as 0.8 fm. The increase in the effective binding 
energy acts to enhance the formation of heavier nuclei. 

The electron fractions {=pe/p^) for the supernova 
matter are compared in the SMB and MFB models as a 
function of baryon density p^ at temperature T =4 and 
8 MeV in Fig. 8. The two results follow nearly the same 
trend (increasingly neutron-rich with increasing baryon 
density); the SMB model, however, predicts somewhat 
less neutron-rich system indicating an effectively lower 
e~-capture rate which may be significant enough to influ- 
ence supernova dynamics [l|, particularly at lower tem- 
perature. In Fig. 9, the calculated baryonic pressure from 
the two models at the aforesaid temperatures are shown 
as a function of density. In contrast to nuclear matter, 
the baryonic pressure for supernova matter in the two 
models are seen to be nearly the same. This is due to the 
fact that in supernova matter, to maintain /3-equilibrium, 
the asymmetry is very high, so much of the matter is in 
the free neutron state that effectively controls the behav- 
ior of the pressure. 

The variation of the symmetry energy coefficients 
(Cf , , and C*°*) of supernova matter as a function of 
density at temperatures T =4 and 8 MeV are shown in 
Fig. 10. These coefficients are evaluated using Eqs. (39) 
and (41) at the relevant asymmetries. The MFB results 
are the same as those of nuclear matter as shown in Fig. 3 
as they are asymmetry independent. In the S'-matrix 
approach to supernova matter, the symmetry energy co- 
efficients have the same trend with density as in nuclear 
matter, however, the values get much reduced. The pres- 



the formation of heavier fragments and hence this reduc- 
tion. 

The baryonic incomprcssibility of supernova matter is 
shown in Fig. 11 as a function of the baryonic density 
in the SMB and MFB models at T =4 and 8 MeV. The 
results in the two models are qualitatively the same at 
both the temperatures and are different from those ob- 
tained for nuclear matter shown in Fig. 5. The calculated 
incompressibility for supernova matter is seen to be sig- 
nificantly larger to that in nuclear matter except at very 
low density where both nuclear and supernova matter are 
mostly composed of free nucleons. 

The baryonic symmetry incompressibilities Kg for su- 
pernova matter in the two models at the aforesaid tem- 
peratures are displayed in Fig. 12. In the mean-field 
model. Kg is independent of asymmetry, hence the MFB 
results are the same as the MF results as shown in Fig. 6. 
In the 5'-matrix approach, the symmetry incompressibil- 
ity for supernova matter at a temperature of 4 MeV is 
much reduced compared to that in nuclear matter at 
lower density; at higher temperature, however, the corre- 
sponding two results are close due to absence of formation 
of heavier fragments within the density range studied. 
The differences between the incompressibilities obtained 
for nuclear matter and supernova matter arise because 
calculations in nuclear matter refer to a fixed asymmetry 
which is taken to be relatively low, but the calculations 
for supernova matter, with changing density, pertain to 
varying asymmetries which turn out to be quite high. 
Furthermore, the presence of the Coulomb contribution 
in the fragment binding energies precludes the formation 
of very heavy fragments, thereby increasing the pressure. 



IV. CONCLUDING REMARKS 

The low-density nuclear or supernova matter is not ho- 
mogeneous, it is clusterized at a finite temperature. In 
this paper, we have investigated some structural proper- 
ties of this low density matter in the the S'-matrix frame- 
work. Special emphasis has been given to understand the 
effects brought out by the isospin structure of the nu- 
clear force. In particular, we have considered its effects 
on the energy, incompressibility and the single-nuclcon 
potentials. The effects are more clearly manifest in the 
symmetry energy coefficients, symmetry incompressibil- 
ity and isovector potential. 

In the S'-matrix approach, symmetric or asymmet- 
ric nuclear matter condenses on isothermal compression 
thereby avoiding the unphysical region brought out in 
the mean-field model. Condensation eases pressure, the 
incompressibility is thus generally lower. Condensation 
has a marked effect on the symmetry incompressibil- 
ity Ks- Around the condensation density. Kg shows a 
peaked structure, particularly at lower temperature, in 
contrast to its smooth monotonic increase with density 
in the mean-field model. In the supernova matter, the 
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are, however, not very different; the rapid neutronization 
in the /3-equiUbrated matter with increasing density is 
mostly responsible for this change. A remarkable feature 
of the symmetry energy in the S-matrix approach for nu- 
clear matter is that as opposed to results from mean-field, 
the symmetry energy is nonlinear in X^. The symmetry 
coefficients are also larger in magnitude in the SM model 
and that the kinetic components are, contrary to expec- 
tations, negative The same trends follow in supernova 
matter. 

In the 5'-matrix framework, the calculated properties 
like the equation of state, the total symmetry energy co- 
efficient and the incompressibilities of nuclear or super- 
nova matter do not depend on any particular choice of the 
nuclear force, they are model-independent. They can be 
directly connected to the experimentally measured phase 
shifts and the binding energies of the fragments consti- 
tuting the matter. The e~-capture envisaged through 
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the chemical equilibrium conditions in supernova matter 
is similarly independent of any model for nuclear interac- 
tions. At low temperatures, in a selective density range, 
the S'-matrix model predicts a somewhat higher electron 
fraction in the system compared to the mean-field cal- 
culations indicating lower e~-capture rate. These may, 
however, play a significant role in supernova dynamics. 



Acknowledgments 

S.K.S. and J.N.D. acknowledge support of DST, Gov- 
ernment of India. M.C. and X.V. partially supported 
by the Consolider Ingenio 2010 Programme CPAN 
CSD2007-00042 and grants FIS2008-01661 from MEG 
and FEDER and 2005SGR-00343 from Generalitat de 
Gatalunya. 



(2009). 

[25] H. Heiselberg and M. Hjorth- Jensen, Phys. Rep. 328, 237 
(2000). 

[26] R. J. Furnstahl, Nucl. Phys. A706, 85 (2002). 

[27] C. Fuchs and H. H. Wolter, Eur. Phys. J. A30, 5 (2006). 

[28] J. Xu, L. W. Chen, B. A. Li, and H. R. Ma, Phys. Rev. 

C 75, 014607 (2007). 
[29] M. Centelles, X. Roca-Maza, X. Viiias, and M. Warda, 

Phys. Rev. Lett. 102, 122502 (2009). 
[30] S. Mallik, J. N. De, S. K. Samaddar, and Sourav Sarkar, 

Phys. Rev. C 77, 032201 (R) (2008). 
[31] J. N. De and S. K. Samaddar, Phys. Rev. C78, 065204 

(2008). 

[32] C. J. Horowitz and A. Schwenk, Nucl. Phys. A 776, 55 

(2006) . 

[33] R. Dashen, S-k. Ma and H.,]. Bernstein, Phys. Rev, 187, 
345 (1969). 

[34] S. L. Shapiro and S. A. Teukolsky, Black holes, White 
dwarfs, and Neutron stars, (John Wiley & Sons. Inc., 
1983). 

[35] A. Bohr and B.R. Mottelson, Nuclear Structure, vol I 

(W.A. Benjamin Inc., 1969). 
[36] R. Dashen and R. Rajaraman, Phys. Rev. D 10, 694 

(1974). 

[37] R. Dashen and R. Rajaraman, Phys. Rev. D 10, 708 
(1974). 

[38] Ronald J. Adler, T. K. Das, and Alvaro Ferraz Filho, 

Phys. Rev C16, 1231 (1977). 
[39] L. I. Schiff, Phys. Rev 133, B802 (1964). 
[40] E. O'Connor, D. Gazit, C. J. Horowitz, A. Schwenk, and 

M. Barnea, Phys. Rev. C 75, 055803 (2007). 
[41] P. Danielewicz, Nucl. Phys. A727, 233 (2003). 
[42] J. N. De and S. K. Samaddar, Phys. Rev. C 76, 044607 

(2007) . 

[43] W. D. Myers and W. J. Swiatecki, Nucl. Phys. A601, 
141 (1996); Lawrence Berkeley Laboratory Report LBL- 
36803. 

[44] C. Ishizuka, A. Ohnishi, and K. Sumiyoshi, Nucl. Phys. 
A723, 517 (2003). 



10 



0.012 



0.008 - 



CO ~0.004|-/ 



> 








(a) 



1 — ' — I — ' — r 
— X=0.0 (SM) 
--X=0.3(SM) _ 



/ 



T=4 



r 



4- 



4- 



X=0.0(MF) 
^ 0.061- - X=0.3(MF) 








0.005 0.01 0.015 0.02 

B -3> 

P (fm ) 



FIG. 1: The isotherms for nuclear matter at T —A MeV 
(upper panel) and T =8 MeV (lower panel) for symmetric 
{X =0.0) and asymmetric (X=0.3) nuclear matter in the S- 
matrix (SM) and mean-field (MF) models. 
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FIG. 2: The single- nucleon {Vn and Vp) and isovector 
{Vn - Vp) potentials in MeV shown at T =4 MeV (left pan- 
els) and T =8 MeV (right panels) for symmetric {X = 0.0) 
and asymmetric {X = 0.3) nuclear matter in the SM and MF 
models as a function of the baryon density. 
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FIG. 3: The interaction, kinetic, and total symmetry energy 
coefficients {Cl,C^ , and Cs°*, respectively, all in MeV) for 
nuclear matter shown as a function of the baryon density at 
T =4 MeV (left panels) and T =8 MeV (right panels) in the 
SM and MF models. 
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FIG. 4: The anharmonicity Dse of the symmetry energy (see 
text) plotted as a function of the baryon density at T =4 and 
8 MeV for three different asymmetries. 
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FIG. 5: The incompressibility coefRcient K shown as a func- 
tion of the baryon density at T =4 and 8 MeV for symmetric 
and asymmetric nuclear matter (left panels) in the SM and 
MF models. In the right panels K is shown as a function of 
temperature in the two models at the baryon densities indi- 
cated. 
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FIG. 6: The symmetry incompressibility coefficient Ks shown 
as a function of the baryon density at T= 4 MeV (upper 
panel) and T= 8 MeV (lower panel) for asymmetric nuclear 
matter with X= 0.1 and 0.3 in the SM and MF models. In 
the MF model, Ks is independent of asymmetry. The ar- 
rows denote the density beyond which the system enters the 
unphysical region in the MF model. 



15 



> 



200 
150 
100 
50 

150 
100 
50 



' 1 
'- (a) 


1 1 1 1 1 

SM(X=0.1) 

SM (X=0.3) 

MF 

p^=0.001 1 


^ \ ^ — =h 

1 (b) 


— ^ — ^ — ^ — 

p''=0.01 

■ 




f - 

1 



8 



T (MeV) 



FIG. 7: The symmetry incompressibility Kg for asymmetric 
nuclear matter shown as a function of temperature at baryon 
densities p-^ =0.001 and 0.01 fm"^ in the SM and MF models. 
In the MF model, Ks does not depend on asymmetry. The 
arrow in panel (b) denotes the temperature below which the 
system enters the unphysical region in the MF model. 
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FIG. 8: The electron fraction Ye in /3-equiIibrated supernova 
matter as a function of the baryon density at T— 4 and 8 
MeV in the S-matrix (SMB) and mean-field (MFB) models. 




FIG. 9: The isotherms for /3-equihbrated supernova matter 
at r=4 and 8 MeV in the SMB and MFB models. 
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FIG. 10: The symmetry energy coefficients Cl,Cf , and Cl°^ 
(in MeV) of /3-equilibrated supernova matter as a function of 
the baryon density for T=A and 8 MeV in the SMB and MFB 
models. 




FIG. 11: The incompressibility of /3-equilibrated supernova 
matter as a function of the baryon density at r=4 and 8 
MeV in the SMB and MFB models. 
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FIG. 12: The symmetry incompressibility Ks for supernova 
matter as a function of the baryon density at r=4 and 8 MeV 
in the SMB and MFB models. 



